function [amplitude, phase] = dft3(x)
    Xk = fft(x); 
    amplitude = abs(Xk);
    phase = angle(Xk);
end